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Within the Standard Model the Higgs couplings to fermions and gauge bosons, as function of the 
particle mass, are predicted to lay on a single line. However, new patterns appear in multi-Higgs 
models that employ the Froggart-Nielsen mechanism, where the diagonal couplings could lay on 
different lines and flavor-violating Higgs couplings could appear too. These aspects are studied for a 
specific model with 3+1 Higgs doublets and a singlet FN field. Constraints on the model are derived 
from the Higgs search at LHC, and some remarks are presented on their implications for the rare 
top and Higgs decay, t ^ ch and h ^ ry. 


I. INTRODUCTION 

The discovery of a Higgs-like particle with ruh = 125 — 126 GeV at the LHC Bia, has verified the mechanism 
of Electro-Weak symmetry breaking of the Standard Model (SM) [3]. Current measurements of the spin, parity, 
and couplings of the Higgs, seem consistent with the SM. On the other hand, several scenarios for Physics 
Beyond the SM (PBSM) have been proposed to address some of its open problems, such as hierarchy, flavor, 
unification, etc HE]. So far, the LHC bounds on its scale (A) are entering into the multi-TeV range, which 
could mean A >> = 246 GeV. A more solid conclusion needs to wait for the next LHG run, with higher 

energy and luminosity. Many papers have been devoted to study the LHC implications for the Higgs couplings, 
for instance in HIT]. The couplings of the Higgs particle to a pair of massive gauge bosons or fermions, are 
proportional to the particle mass. However, the LHC has tested only a few of these couplings, i.e. the ones with 
the heaviest SM fermions and W, Z. Non-standard Higgs couplings, including the flavor violating (FV) ones, 
are predicted in many models of physics beyond the SM, for instance in the general multi-Higgs models 0191 
and SUSY |33]. 

Many ideas have been proposed to adress the flavor problem m, for instance: Textures and GUT-inspired 
relations, flavor symmetries and radiative generation. Within the flavor symmetry approach, the Froggart- 
Nielsen mechanism assumes that above some scale Mp, such symmetry forbids the appearence of the Yukawa 
couplings; SM fermions are charged under this symmetry (which could be of Abelian type However, 

the Yukawa matrices can arise through non-renormalizable operators. The Higgs spectrum of these models 
could include a light and heavy Higgs boson. In these models the diagonal flavor conserving (FC) couplings of 
the light SM-like Higgs boson could deviate from SM, while flavor violating (FV) couplings could be induced 
at small rates too. Within the SM, the FG Higgs couplings to fermions and gauge bosons, as function of the 
particle mass, lay on a single straight line. However, in multi-Higgs models, they could lay on distinct lines. As 
illustrated in Figure 1, the scalar spectrum could include: 

i) A lightest state that should be identified with the SM-like Higgs boson. One expects, the appearence of 
small deviations from SM predictions for the diagonal Higgs-fermion couplings and suppressed FV couplings. 

ii) States with Flavon-dominated composition, which could provide the more radical signature of the models 
under consideration. The observation of these signals depends on the flavon scale, and could at the reach of the 
LHG reach if such scale were about 0(1) TeV. 

iii) Heavy Higgs bosons which could have large mixing with flavons, and thus deviate significantly from SM 
expectations for FV couplings, that could also be searched at LHG. 


II. HIGGS COUPLINGS WITHIN A 3+1 HIGGS MODEL 

Thus, we shall consider a 3+1-Higgs doublet model, denoted as fl>o, 5 ^ 3 - The fl>i gives masses to 

up-type quarks, while fl >2 and fl >3 give masses to d-type quarks and leptons, respectively. We shall impose a 
discrete symmetry in such a way that one doublet (fl>o) is of the inert-type, and therefore contains a Dark matter 
condidate [12]. Furthermore, we shall also include Froggart-Nielsen scalar field (SM singlet S). The possibility 
of having light flavon fields was studied in ref. m, and more recently in The Yukawa lagrangian is 
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FIG. 1: The Higgs and flavon spectrum in multi-Higgs models. 


given by: 


+ h.C. 


( 1 ) 


where n,p,q denote the charges of each fermion type under an Abelian flavor symmetry, which will help to 
explain the fermion mass hierarchy. The flavon field S is assumed to have flavor charge equal to - 1 , such that 
Cy is invariant. Then, Yukawa couplings arise after the spontaneous breaking of the flavor symmetry, 

i.e. Aa; = where < S > denotes the flavon vacuum expectation value, while Mp denotes the heavy 

mass scale, wfiich represents the mass of heavy fields that transmit such symmetry breaking to the quarks and 
leptons. 

The Higgs and Flavon fields are written in terms of mass eigenstates, through the rotation (4x4): 

ReS = ( 2 ) 

Furthermore, as the vev’s must satisfy: = v‘^, with v = 246 GeV, we find convenient to use spherical 

coordinates to express each vev (vi) in terms of the total vev v and the angles jli and (32, as shown in figure 2 , 
namely: Vi = v cos Pi, V 2 = v sin Pi cos P 2 and Vs = v sin Pi sin P 2 . 





FIG. 2: The Higgs vevs in spherical coordinates. 

Thus, for the lightest Higgs state {hi = h^, one gets finally the following interaction lagrangian for the 
Higgs-fermion couplings, 

Cy = [—UM^U + l^DMdD + 1-LMiL + nWiZ^Uj + K'^DiZ'^D. + K^iZnAh^ (3) 

V V V 

where: 77 ^ = Of^/cos/^i, 77 ^ = O 21 / A cos ^^ 2 , r]^ = / sin sin ^^ 2 , describe the strength of the flavor- 

diagonal Higgs couplings. While the FV Higgs couplings are described by the parameters: rN = -Oji cos/3i, 
= ^Oji sin Pi coS;d 2 , cos^di sin;d 2 . 
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Besides the Yukawa couplings, we also need to specify the Higgs couplings with vector bosons, which is written 
as ghvv = XvOhVV^ factor xv given as: 


XV - —On + —O21 + —O31 

= cos /3i Oji + sin /3i cos /32 O^ + sin /3i sin ^2 Oji 


( 4 ) 


It is interesting to note that the coupling xv can be written in terms of the FC fermionic couplings, which 
can be seen as a type of sum rule, i.e. 


yy = cos^ /3i + sin^ pi cos^ 132 T + sin^ pi siiP P 2 rf 


( 5 ) 


Moreover, since the Higgs couplings to first generation quarks and leptons is highly suppressed, in order to 
study the FV Higgs coupling, which depends on the matrices , we shall consider the 2-3 family sub-system. 
Namely, for up quarks the Z-matrix (in mass eigenstate basis), is given by: 




\ru \^u \ 
^22 -^23 \ 

4?3 ) 


(6) 


and similarly for d-quarks and leptons. We find a relation among the parameters, such that we can express 
the in terms of the ratios of masses and the CKM angle Vcb — 523- Namely, we define: ru = niclmt, rd = 

rris/mb, and = ¥^ 2 /^ 33 ^ ^’2 = ^ 23 / 433 - Similarly: rf = ¥^ 21 ^^ 33 , = ¥^ 3 /¥^ 3 . Within this approximation 

we have: Fgg ~ F 33 for f = u,d. Then, r( = rf + r^, and the ratios of Yukawas must satisfy the following 
relation: 


ri = 

1 + 


<523 


( 7 ) 


Thus, in order to study the predictions of our model, we need to specify the vevs Vi and the rotation matrix 
for Higgs particles (Oij). For the 2 HDM (see for instance [IS]), LHC Higgs data favors both decoupling and 
alignment solutions, namely both tan^d 1 and tan^d >> 1 are acceptable solutions. Thus, for the vev’s we 
leave pi as free parameter, then explore the following cases: 

• (VEVl) We can take first V 2 = V 3 , which in spherical coordinates, means: P 2 = 

• (VEV 2 ) We also consider vevs with V 2 < '^ 3 , for which we take: P 2 = 

• (VEV 2 ) We also consider vevs with V 2 > '^ 3 , for which we take: P 2 = 

Then, for the rotation matrix O of real components of scalar fields, we can identify several interesting scenario 
where the 126-Higgs is lighter than the heavy Higgs particles and the fiavons, i.e. ruh < m//. nfiHF^ which 
have masses of order TeV. Here, we shall consider a special sub-case, namely we shall assume that > O^, 
and will use the orthogonality relation for the rotation matrix O, in order to relate the parameters, namely: 

+ (OH)^ + {O 31 Y + = 1- Eurthermore, assuming Oj^ (for i ^ j) one has Oj^ = 


III. FLAVOR CONSERVING HIGGS COUPLINGS AT LHC 

LHC data on Higgs boson has been used to derive bounds on the Higgs couplings, i.e. deviation from the 
SM, which are defined as: gnxx = Ohxxi^ 3- ex)^ where 77 ^ = 1 -h ex- A complete analysis is done in ref. [7]; 
for fermions, they obtain the allowed values: et = — 0.21 ± 0.23, e^ = —0.19 ± 0.3, e^- = 0 ± 0.18; while for W 
(Z) bosons they find: ew = —0.15 ± 0.14, ez = —0.01 ± 0.13. 

An extensive analysis of parameters satisfying these bounds will be presented elsewhere, with detailed numer¬ 
ical scans; here we shall pick a few specific points in parameter space, which satisfy the LHC bounds, and will 
help us to understand qualitatively the beahivour of the model. These poins will also be used in the next section 
in our analys of FCNC top decays. Thus, we show in figure 3-6 the predictions for each of these parameters, 
as function of the angle Pi, for the case with ^^2 = |, f, f, and for On = 0.5, 0.75, 0.9. We can see that it is 
possible to satisfy these bounds for all the e’s. 
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FIG. 3: Higgs coupling factors et and ez for the parameters of set 1. The horizontal lines are the experimental limits 
on these factors. 
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FIG. 4: Higgs coupling factors Cr and for the parameters of set 1. The horizontal lines are the experimental limits 
on these factors. 


One specific point, in agreement with all data, is : Pi = 0.5 with On = 0.9 and /32 = |. For these values 
we have: rf = 1.03, 77 ^ = 0.6 and rf = 1.04 and Xv = 0.96. This shows that h behaves very much SM-like, 
except for the coupling with d-type quarks. Then, using these values we can plot the Higgs-fermion coupling as 
function of the mass, as shown in figure 2. We can see that the couplings for each fermion type lay on different 
lines, which could be distinguished from the SM (Black line). Future meassurements of these couplings at LHC 
Run2, ILC or FCC will help us in order to discriminate betwen our model and the SM. 



FIG. 5: The Higgs-fermion coupling factors as function of mass, for parameters defined in the text. SM case (Black), 
up-type quarks (red), d-type quarks (blue), charged leptons (green). 


We also find that the corrections contained in the factors , will not change significantly the above dis¬ 
cussion for the top quark-Higgs couplings. However, the Higgs coupling with the lighter fermions (56, cc, r+r“), 
could be measured at next-linear collider (NLC) with a precision of a few percent, and it will be possible to test 
these effects. The corrections to the coupling hbb, could modify the dominant decay of the light Higgs, as well 
as the associated production of the Higgs with b-quark pairs [16]. 

We shall consider the following sample values: = 0.05,0.1,0.3, and also assume: cos/3i 1, then table 

1 shows the values of the entries for the Z'^ matrix for the 2nd-3rd family case. We choose to focus on the 
up-quark sector, because we want to get an estimate for the most relevant predictions of the model, which we 
believe is related with the top quark physics, and in particular for the decay t cPh. For the specific point 
in parameter space, presented in previous section: Pi = 0.5 with On = 0.9 and P 2 = which is in agreement 
with LHC data, we obtain the following value rN = 0.23^. 
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Scenario 

u[TeV] 

X Z 23 

BR.{t ch) 

XI 

0.5 

1.2 X 10“^ 

8.6 X 10“® 

X2 

1 

6.1 X 10“® 

2.2 X 10“® 

X3 

10 

6.1 X 10“® 

2.2 X 10““ 

Y1 

0.5 

6.9 X IQ-'^ 

2.7 X IQ-'^ 

Y2 

1 

3.4 X IQ-'^ 

6.8 X 10“® 

Y3 

10 

3.4 X 10“^ 

6.8 X 10“® 

Z1 

0.5 

2.9 X 10“^ 

4.8 X 10“'‘ 

Z2 

1 

1.4 X 10“2 

1.2 X 10“'‘ 

Z3 

10 

1.4 X IQ-'^ 

1.2 X 10“® 


TABLE I: The factor x Z 23 and Branching ratios for t ^ ch 


IV. THE FCNC DECAY t ch 

The top radiative decay t ^ c + 7 was first calculated in m, followed by ref. m which presented the 
complete calculations of the FCNC modes t cX {X = 7 , g, Z, h); for corrections to SM results for t ^ ch and 
SUSY results see also: [191 [20]. The 3-body FCNC decays modes t cW^W~{Z were presented in ref. 
[2T1I22], while the mode t cl ~was discussed recently in [23] [24]. The 4-body decay t bWl~t^ was also 
studied recently m- Here we shall focus on the mode t ch, which can reach large BR’s. The decay with for 
t ^ ch within our model is given by: 


r(i^c/l) = ^|K“Z23p (8) 

Using the value r(t ^ 6 + W) ^1.5 GeV, we obtain: BR{t ch) = 0.58|/^^Z23p. For vju = 0.25, and one 
finds that the B.R. could reach a value BR ^ 1.5 x 10“^, which could be tested at LHC [26|. Values of BR for 
other choices of parameters are shown in table 1 . 


V. THE LEV DECAY h^ry 

Another interesting probe of FV Higgs couplings is provided by the decay h ry, which was initially studied 
in refs. [271128]. Subsequent studies on detectability of the signal appeared soon after [29]-[3T]. Precise loop 
calculations with massive neutrinos, SUSY and other models appeared in [321435] . The recent search for this 
decay at LHC m, have resulted in a bound for the corresponding branching ratio of order B.r.{h ry) < 
1.51 X 10“^ at 95% C.I.. Furthermore, given that the best fit to the data gives B.r.{h ry) = O.SdlQ'Jy x 10“^, 
many more papers have appeared recently, trying to explain this result [38]. The search for this LFV Higgs 
decay could be one great opportunity to find new physics at the LHC Runll. 


VI. CONCLUSIONS AND OUTLOOK 

We have studied the Higgs couplings, within a model with 3+1 Higgs doublets, where the masses for each 
fermion type, arise from a different Higgs doublet. This model also includes mixing of the Higgs doublets with a 
Flavon field, which generates the Yukawa hierarchies and induces Flavor-violating Higgs couplings at acceptable 
rates. Constraints on these couplings, derived from Higgs search at LHC, and their implications for FCNC top 
decay t ch, were discussed too. We find that this mode could reach a BR of order 10“^, which could be 
studied at LHC. In the down-quark and lepton sectors, there are also interesting aspects to study in the future, 
such as the rates for rare b-decays. or the decay h ^ ry, which can be induced at rates that could be detected 
at future colliders. 
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